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Abstract 

The pressure in a classical Coulomb fluid at equilibrium is obtained from 
the Maxwell tensor at some point inside the fluid, by a suitable statistical 
average. For fluids in an Euclidean space, this is a fresh look on known 
results. But, for fluids in a curved space, a case which is of some interest, 
the unambiguous results from the Maxwell tensor approach have not been 
obtained by other methods. 
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1. INTRODUCTION 



The aim of the present paper is to revisit the concept of pressure for a Coulomb 
fluid, i.e. a fluid made of particles interacting through Coulomb's law (electrolyte, 
plasma,...). We consider only fluids in thermodynamic equilibrium, and assume that 
classical (i.e. non-quantum) statistical mechanics is applicable. 

Pressure is often defined as the force per unit area that a fluid exerts on the walls 
of a (large) vessel containing it. However, pressure may also be defined without 
reference to any wall. One has to imagine some immaterial plane surface across the 
fluid, and pressure is the force per unit area with which the fluid lying on one side 
of this surface pushes on the fluid lying on the other side. Both definitions agree 
with each other. 

Prom a microscopic point of view, the force between two parts of the fluid is 
usually described in terms of the interactions between the molecules. Two molecules 
at a distance r from each other are supposed to interact through some potential. 
This is the standard approach, which is briefly recalled in Section 2. 

In the case of electromagnetism. Maxwell, following Paraday, came to a different 
point of view: the forces between two charged objects are mediated by fields. At 
a given point of space, even in vacuum outside the charges, there is a stress tensor 
(the Maxwell stress tensor), a local quantity defined in terms of the fields at that 
point, similar to the stress tensor within some elastic medium. In this picture, every 
region of "empty" space exerts forces on the regions ajacent to it. In Section 3, it 
is shown how the pressure at some point inside a Coulomb fluid can be defined and 
computed from the Maxwell tensor at that point by a suitable statistical average, 
with an appropriate prescription for obtaining a finite result. This result agrees with 
the standard one. 

Section 3 "extends" the above ideas to the case of two-dimensional models. 
Section 4 discusses the case of Coulomb fiuids living in a curved space. In this 
case, the Maxwell tensor approach will be shown to be especially appropriate. 
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2. A SUMMARY OF THE STANDARD APPROACH 



In the simple case of a fluid made of one species of particles, with a pair inter- 
action f (r) depending only on the distance r, the pressure P is found to be, in the 
thermodynamic limit. 



where n is the number density (number of particles per unit volume), k is Boltz- 
mann's constant, T is the temperature, and g{r) is the pair distribution function. In 
(2.1), nkT is the ideal gas part of the pressure (related to the momentum carried by 
the particles), while the following term, due to the interactions, is called the excess 
pressure P^x- The same equation (2.1) is obtained by looking either at the pressure 
on the walls*^^^ or at the pressure in the bulk fluid*^^-*. 

Real Coulomb fluids are made of several species of particles (for instance, in 
an electrolyte, positive and negative ions, plus the solvent molecules) . In a classical 
model, some short-range non-Coulombic interaction must be introduced, for avoiding 
the collapse on each other of particles of opposite sign. Here, for simplicity, we rather 
consider only a simplified model, the one-component plasma (OCP)^^-*: identical 
point-particles of one sign, each of them carrying an electric charge g, embedded in 
a uniform background of opposite sign which ensures overall neutrality. Only the 
Coulomb interaction is retained, thus the interaction is v{r) = q'^/r and {dv/dr)r = 
—q^/r. Due to the background, the average charge in a volume element dr at a 
distance r of a given particle is qn[g{r) — 1] rather than qng{r), and, in the case of 
the OCP, equation (1) is to be replaced by 



where h{r) = g{r) — 1 is the pair correlation function. It may be noted that Pex is 

one third of the potential energy density.[] 

The above standard approach is based on the assumption of an interaction-at- 

distance g^/r. In the next Section, it will be shown how (2.2) can be derived by 

using the Maxwell tensor. 

^For an OCP, there are several non-equivalent possible definitions of the pressure. The pressure 
(2.2) is the thermal pressure, in the sense of Choquard et al.^'*-' 
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3. THE MAXWELL TENSOR APPROACH 



If only electrostatic interactions are retained (magnetic effects are neglected), 
the Maxwell tensor is^^^ 

Tafi = -T-{EaEl3 - • E5a/3) (3.1) 

In (3.1), the Greek indices label the three Cartesian axes {x,y,z). T^p is the a 
component of the force per unit area transmitted, across a plane normal to the j3 

axis, to the fluid lying on the negative side of this plane. Thus, choosing for j3 any 
axis, say the x axis, one obtains for the excess pressure, which is a force along that 
axis, 

Pe. = - < >= < - - > (3.2) 

where < . . . > denotes a statistical average on all particle configurations (the electric 
field at some point is a function of the particle configuration) . Our task is to evaluate 
the statistical average (3.2) at some point inside the fluid, say at the origin. 

Let p*-^-'(ri2) be the statistical average of the product microscopic charge density 
at ri times microscopic charge density at r2 (ri2 = r2 — ri). From (3.2), 

In the present case of an OCP, 

P^^\ri2) = g'H(ri2) + n^h{n2)] (3.4) 

Using (3.4) in (3.3) gives to P^x two contributions Pgeif and Pnonseif involving the 
5 part and the h part of (3.4), respectively. Pnonseif gives no difficulty. This is 

a convergent integral (indeed, /i is -1 at small because the particles strongly 
repel each other, and h has a fast decay at large because remote particles are 
uncorrelated) . Because of the rotational symmetry around the origin, it can be 
rewritten as 

Pnonseif = |^ / dridr2^^^j^h{ru) (3.5) 

But 

Stt y r6 ^ ' 
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self 



diverges at small r. 

The resolution of the difficulty is that the force that each particle exerts on itself 
should not be taken into account. Thus, the integral in (3.6) must be regularized by 
the prescription that no particle sits on the x — plane on which we have chosen 
to compute the pressure force. This prescription can be enforced by removing from 
the integration domains a thin slab —e<x<e and taking the limit £ ^ at the 
end. This prescription does not change the convergent integral (3.5). But it means 
that the self part (3.6) must be defined, in cylindrical coordinates {x,p), as 



Pseif ^ -^^i^^l dx 27Tdpp (3.7) 
Stt £^oJ|2;|>£ Jo [x^ + p^y 

Since the integral on p, performed first, is found to vanish, the result is Pgeif — 0. 

As to Pnonseif, (3.5) Can be easily computed by taking as integration variables 
ri and ri2, and performing the integral on ri ffist with the result 47r/ri2. The final 
result is 

r 1 



Pex = Pnonseif = / h{ri2) (3.8) 

6 J ri2 
in agreement with the standard formula (2.2). 

An alternative way of calculating Pseif will turn out to be more appropriate for 
extensions which follow. (3.6) is split into the contributions Pq oi r < ro and Pi 
of r > ro, where ro is the radius of a small sphere centered at the origin. The 
prescription that no particle sits on the plane x = does not change the convergent 
part Pi, which can be computed, using the rotational symmetry, as 

= ^ = ^ (3.9) 

247r Jr>ro 6ro ^ ^ 

It is only in Pq that the rotational symmetry is broken by the prescription \x\ > e, 
which gives 



"''lim/ ,,f^2ndpp4^,^-'^ (3.10) 



Stt £^oJs<\x\<ro Jo (a;2 + p2)3 6ro 

Thus Pseif = Po + -Pi = 0, and (2.2) is retrieved. 
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4. TWO-DIMENSIONAL MODELS 



Two-dimensional models of Coulomb fluids are of interest for at least two rea- 
sons. First, some of these models are physically relevant. Second, exact results are 
available. The two-dimensional case has special features which require the present 
separate discussion. 

In two dimensions, the Coulomb interaction (as defined through the Poisson 
equation) between two charges q and q' is —qq'\n{r/L), where L is some irrelevant 
length. Since this interaction diverges at r = only mildly, in addition to the OCP 
it is also possible to consider a two-component plasma (TCP), made of positive and 
negative point-particles of respective charges q and —q, without any additional short- 
range repulsion (which is stable provided that the coupling constant F = q"^/ kT be 
smaller than 2). 

For the OCP, the two-dimensional analog of (2.1), with the background taken 
into account, is 

P^nkT- -n^ [ ^rh(r)dr (4.1) 
A J dr 

Now v{r) = —q^hi{r/L) and (4.1) becomes 

P = nkT + ^n^q^ J h{r)dv (4.2) 

Perfect screening, present in a conductor, says that 

n j h{r)dv = -1 (4.3) 

(this means that the polarization cloud around a particle of charge q carries the 
opposite charge —q). Using (4.3) in (4.2) gives the simple exact equation of state ^^'^^ 

P = n{kT - ^) (4.4) 

Now, Pex = —nq'^/4 is no longer related to the potential energy density. 

We now turn to the Maxwell tensor approach. In two dimensions, the Maxwell 
tensor is 

Tap = ^{EaEp - ^E • E5af3) (4.5) 
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with Greek indices now labeling two Cartesian axes {x,y). (3.3) is replaced by 

Pex = --r / "ridrs '{ri2) (4.6) 

47r J ri r2 

where, for an OCP, (3.4) still holds. Now, although (4.6) still converges for large 
values of ri and r2 (because p^'^\r 12) has a fast decay as ri2 increases and its integral 
vanishes), separating it in self and nonself parts would generate terms separately 
diverging at infinity. Here, it is more appropriate to split (4.6) in another way, 
similar to what has been done at the end of Section 3. Namely, in (4.6), one 
separates the contribution Pq of the integration domain (ri,r2 < Tq) and the rest 
Pgx — Po- This rest is a convergent integral and, by rotational symmetry, it vanishes. 
One is left with Pq which can be split into its self and nonself parts, with now a 
nonself part which is convergent and also vanishes by rotational symmetry. Finally, 
the self part has to be defined in the same way as (3.10), and 



Thus 



Pq — — : — lim / dx ^ dy—- = — (4.7) 

47r ^-oA<|a:|<ro J-y/^J^ ^ {x^ + y^f 4 ^ ^ 



Pex = Po = (4.8) 



in agreement with (4.4). 

Similar considerations hold for the TCP, as long as F < 2, and the equation of 

state again is (4.4), where now n is the total number density of the particles. 

The two-dimensional OCP can also be obtained as a limit of the i/- dimensional 
one, as explained in Appendix A. 



5.CURVED SPACES 



The statistical mechanics of a Coulomb fluid living in a curved space is of interest 
for at least two reasons. First, for doing numerical simulations (necessarily on a finite 
system) without having to deal with boundary effects, a clever method has been to 
confine the system on the surface of a sphere (in the two-dimensional case)^^^ or 
an hypersphere (in the three-dimensional case)*^^'^°\ Second, for two-dimensional 
Coulomb fluids on a surface of constant negative curvature (pseudosphere) it 
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is possible to go to the limit of an infinite system while keeping a finite curvature, 
thus to look at the properties of a curved infinite system (something which cannot 
be done for a sphere or hypersphere) . 

The present paper actually arose from the question: How to define the pressure 
of a Coulomb fluid in a curved space, away frow any wall? A formula like (2.1) is 
based on the interaction-at-distance picture: the force acting on the fluid lying on 
one side of some immaterial plane is the sum of elementary forces acting on each 
molecule. This picture cannot be generalized to the case of a curved space, because 
there is no straightforward way of summing forces (vectors) applied at different 
points of space. Thus, the Maxwell tensor picture seems to be the only possible one, 
defining the pressure at a given point of space as a local quantity depending only 
on the electric field at this point. 

Three kinds of Coulomb fluids will be considered: The three-dimensional OCP on 
a hypersphere, the two-dimensional OCP or TCP on a sphere, the two-dimensional 
OCP or TCP on a pseudosphere. 

5.1. OCP on a Hypersphere 

The hypersphere is the four-dimensional analog of the usual sphere. We consider 
an OCP living on the three-dimensional "surface" 6*3 of a hypersphere of radius 
R. On 5*3, the geodesic distance between two points is Rip, where G [0, vr] is the 
angular distance between these points, as seen from the center of the hypersphere. 
The volume element between two concentric spheres of radii Rip and R{ip + dip) is 
dV = AttR^ sin^ ipdip. The total volume of S3 is V ^ 2t:'^R^. 

Since S3 is a compact manifold without boundary, electric potentials and fields 
can be defined only if the total charge is zero. In particular, the electric field 
created by one point charge cannot be defined. For overcoming this difficulty, one 
can consider the OCP as a collection of pseudocharges '^^^: a pseudocharge is defined 
as a point charge q plus a uniform background of total charge —q. At a point M 
located at a geodesic distance Rip from a pseudocharge located at Mq, the electric 
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potential created by the pseudocharge is 



$= ^ ('(7r--0)ctnV'--) (5.1) 



nRV ' 2 
where Vq is an arbitrary constant. The corresponding electric field at M is 

E= (ctn^ + 1 t (5.2) 

where t is the unit vector tangent to the geodesic MMq at M. Prom (5.1), the 
interaction energy between two pseudocharges i and j at a geodesic distance it^^'ij 
of each other is found to be 

" TTfl (^^ ~ V'uOctn - ^ (5.3) 

independent of Vq. 

The excess pressure is given by (3.2) where the electric field can be written as 
E = X] Ej, with Ej the field created by the ^-th pseudocharge. As above, (3.2) can 
be split into a self part Pgeij (made of EjEj terms) and a nonself part Pnonseif (made 
of EjEj- {i^ j) terms). 

Because of the rotational symmetry, Pnonsei f can be written as 
„ 1 1 

Ei ■ Ej) = -Unonself (5.4) 



nonself 



247r^^ ' " 3 



where Unonseif is the nonself part of the potential energy density, which can be 
reexpressed in terms of the interaction rather than in terms of fields, by the usual 
integration by parts, as 

Unonself = ^-JBjdl^ (^ij)) = y / <PWh{lp)dV (5.5) 
i<j 

where one can use the pair correlation function h rather than the pair distribution 
function g — h + 1 since / (f)dV — 0. 

As to Pseif, it is a divergent integral which however can be made finite by adapt- 
ing what has been done at the end of Section 3, namely splitting it into the contribu- 
tion Pq of geodesic distances Rip < Ripo, and the finite contribution Pi of Rip > Ripo 
for which rotational symmetry can be used. In terms of E{ip) given by (5.2), 

Pi = — / E2(^)4rfsin^'0# (5.6) 
2471 Jvo 



Evaluating the integral in (5.6) as V'o — >^ gives 

On the other hand, for Pq, the curvature effects become negligible as ipo ^ (more 
precisely, as shown in Appendix B, they are 0{ijjo)) and the Euclidean prescription 
(3.10) can be used, with tq = Ripo, giving 

^° = -6^ + «(^°) (^■«* 
Finally, the total pressure is 

2 2 

P^nkT+'^l <P{^|;)h{^|;)dV - ^ (5.9) 

This is the generalization of (2.2) to the case of a hypersphere. 

The present evaluation of the pressure makes no exphcit use of the self-energy of 
a pseudoparticle. This is an important remark. Indeed, another possible definition 
of the pressure would be minus the derivative of the free energy with respect to the 
volume. But, for evaluating the free energy, it is necessary to define properly the 
zero of energy for a system of pseudoparticles, and this necessarily involves some 
heuristic convention about the self-energy of a pseudoparticle. In ref.9. a reasonable 
convention gave —Sg'^/AnR for this self-energy, and the corresponding pressure is 
identical with (5.9). This pressure does obey the usual relation 

Pex = lu (5.10) 

where u is the total potential energy density defined with the above convention. 
However, in ref.lO, another reasonable convention (giving a faster approach to the 
thermodynamic limit as i? ^ oo) has been used, with additional terms of order 
higher than 1/R, and the pressure derived from the corresponding free energy no 
longer agrees with (5.9). The definition of the pressure in terms of the Maxwell 
tensor is free from this arbitrariness. 
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5.2.0CP or TCP on a Sphere 



The above considerations can be easily adapted to the (simpler) case of two- 
dimensional Coulomb systems living on the surface of a sphere. Now, the electric 
potential created by a pseudocharge is 

$ = -gin sin (^/2) + 1/0 (5.11) 

Pnonseif ^-nd Pi, Convergent integrals, vanish because of the rotational symmetry. 
One is left with Pq, for which the curvature effects are negligible and (4.7) holds 
with the same result (4.8) as in the case of a plane system. This result Pex — —nq^/A 
holds for an OCP, and also for a TCP when F < 2. 

Here too, the free energy depends on an arbitrary convention about the zero of 
energy. In ref.8, this convention was implicitly made by the way in which (5.11) was 
used together with the choice Vq = —qln{2R/L), with R the radius of the sphere. It 
is only thanks to this convention that the corresponding free energy has a derivative 
with respect to the sphere area which correctly gives the equation of state (4.4). 



5.3.0CP or TCP on a Pseudosphere 



Recently, two-dimensional Coulomb systems living on a surface of constant neg- 
ative curvature (a pseudosphere) were studied. ^^^-^ Since, unlike a sphere, a pseudo- 
sphere infinite, one has the interesting possibility of considering systems which are 
both infinite and curved. 

Let a be a length such that the Gaussian curvature of the pseudosphere is — 1/a^ 
(instead of 1/i?^ on a sphere). Now, the electric potential and field created by a 
single point charge q exist. At a geodesic distance s from this charge, the electric 
potential is 

$ = -glntanh^ (5.12) 

where the possible additive constant has been fixed by the condition that this po- 
tential vanishes at infinity (s — > 00). The electric field is 

E - — TTTTt (5.13) 
asmh(s/aj 
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where t is the unit vector tangent to the geodesic. 

The pressure can be obtained from the Maxwell tensor just as in the case of a 
sphere, with the same result (4.4), i.e. Pex — —nq^/A. This pressure holds not only 
in the thermodynamic limit, but also at the center of a finite disk. 

The above result for the pressure calls for some discussion. On a pseudosphere, 
when the size of a large domain increases, its perimeter grows as fast as its area. As 
a consequence, there is no uniquely defined thermodynamic limit for the free energy 
per particle (this limit depends on the shape of the domain and on the boundary 
conditions). A bulk pressure cannot be defined by deriving the free energy with 
respect to the area. In rcf.ll, a bulk "pressure" p was defined by its virial expansion 
with the prescription that the thermodynamic limit of each virial coefficient Bk 
(which seems to exist on a pseudosphere) has to be computed before the virial series 
in powers of the density n is summed. It is now apparent that this p is not identical 
with the pressure P obtained from the Maxwell tensor in the form (4.4). We now 
believe that the correct pressure is P, while p only is a mathematical quantity 
(seemingly with some interesting properties). 

Nevertheless an important result of ref.ll is true: there is at least one thermody- 
namic quantity, the bulk energy per particle, which has a scries expansion in integer 
powers of the density, in contrast to the case of a plane system in which the energy 
per particle is singular at zero density. 

5. 4. Why no Trace Anomaly? 

Some time ago, it has been remarked that conducting Coulomb systems are 
critical-like at any temperature^^^'^^^ in some sense: they have long-range electric 
potential and field correlations, and the free energy of a two-dimensional Coulomb 
system with an electric potential (f) has logarithmic finite-size corrections similar 
to the ones which occur*^^^) in a critical system described by a conformal-invariant 
field theory with a field 0. In ref.l4, it was shown that, for a critical system, this 
logarithmic correction to the free energy is associated to a trace anomaly of the stress 
tensor, proportional to the Gaussian curvature of the surface on which the system 
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lives. Thus, at first sight, one might expect that the pressure of a two-dimensional 
Coulomb system (which is minus one half of the expectation value of the trace of the 
Maxwell tensor, with a suitably defined self part) would have a term 0{1/E?) on 
a sphere and 0(l/a^) on a pseudosphere. Yet, such terms are not present in (4.4). 
Why? 

Actually, for a Coulomb system made of N particles on a sphere of radius R, thus 
with a density n — N/At^B?, the free energy F has a term [kT /&)\n.N . The pressure 
is (with a suitable definition of the zero of energy) a partial derivative at constant 
N: P = — ((?F/(5(47ri?^))jy. Thus the In term in F gives no contribution to the 
pressure, in agreement with (4.4). However, in field theory, some ultraviolet cutoff (a 
length ri) has to be introduced and the trace of the stress tensor has an expectation 
value < © > related to {dF/d{i.T:B?'))^ with now a derivative taken at constant 
cutoff. When the Coulomb system is described in terms of a field theory, the role of 
the cutoff is played by the microscopic scale 77 = n"^/^. Thus, the InN — ln(n47ri?^) 
term in F is associated to a trace anomaly in the field-theoretical < © >, not in P. 

A related statement is: If the expectation value of the trace of the Maxwell 
tensor is computed with a field-theoretical measure (the functional integral measure 
T>(f)), it has a trace anomaly. This trace anomaly is not present when the measure 
is the particle configuration space one dYidY2 . . .drj^. 

Similarly, the pressure (4.8) on a pseudosphere has no trace anomaly. 

6.CONCLUSION 

The pressure in a Coulomb fiuid has been defined as minus the statistical average 
of a diagonal element, say T^x of the Maxwell tensor. This definition leads to an 
ill-defined integral, which however can be given a definite value by an appropriate 
prescription: the fiuid is supposed split into two regions separated by a thin empty 
slab normal to the x-axis, T^x is computed at a point inside this slab, and the limit 
of a slab of zero thickness is taken at the end. 

For Coulomb fluids in an Euclidean space, this approach through the Maxwell 
tensor is just a fresh look on well-known results. But, for Coulomb systems in a 
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curved space, we are not aware of any other way of obtaining an unambiguous value 
for the pressure. 

For simphcity, only point-particle systems without short-range forces have been 
considered. But an extension to systems with hard cores seems feasible. 



APPENDIX A. THE 1/-DIMENSIONAL OCP 



In this Appendix, the excess pressure of a i^-dimensional OCP {u > 2) is related 
to its potential energy density. The dimension i> is treated as a continuous variable, 
and the limit u 2 is taken. 

The unit of charge is defined such that the electric field at a distance r from a 
unit charge be l/r^~^. Thus, the potential is — 2)r''~^. The Maxwell tensor is 

Ta0 = ^{E^Ep - ^E • E (A.l) 

where 

is the area of the sphere of unit radius. 

In terms of the Maxwell tensor T, the nonself part of the pressure is 

Pnonself <i tlC T ^nonself TTc ^ ^nonself (-^•'^) 

where the nonself part of the electrostatic energy density is 
1 ^^2^ <fn^ f , , ri ■ r2 



< E^ >nanself= / dT,dr2^—^h{ri2) (A.4) 



Taking as integration variables ri and ri2, and performing first the integral on ri, 
one finds, as expected, the potential energy density 



^ -Ej -^nonself 



2b,y-i 2 J {i^-2)ri2 

As to the self part of the pressure Pseif — — < T^x >seif, it must be defined, like 
in (3.7), as 
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Here too, the integral on p, performed first, is found to vanisli, tlius Pgeif = 0. 
Tlierefore, tlie final result for the excess pressure is 

g^n^ f 1 



2u ^ ,12 



/ c?ri2^;^/i(ri2) (A.7) 



In the limit i/ — > 2, using the perfect screening rule (4.3), one retrieves 

Pe.--^ (A.8) 



APPENDIX B. CURVATURE EFFECTS IN A SMALL SPHERE 

In this Appendix, (5.8) is derived. 

In four-dimensional Euclidean space, with Cartesian coordinates {x,y,z,t), the 
surface 5*3 of a hypersphere of radius R centered at the origin usually is parametrized 
by the hyperspherical coordinates {u,v,w) related to the Cartesian ones by 

a; = sinwsinwcosM, y = sinw sinv sinu, z = smwcosv, t = cosw 

0<'u<27r, 0<'t;<7r, 0<'u;<7r (B.l) 

However, here, it is more convenient to parametrize S3 by the three independent 
variables {x,y,z). We define r — {x'^ + y'^ + 2;^)^/^ and p — {y'^ + 2;^)^/^. A useful 
relation is = R"^ sin^ w. Using the Jacobian for the change of coordinates, one 
finds that the volume element dV — R^ sin^ w sin vdudvdw becomes 

dV = (B.2) 
cosw 

The hypersphere pole x = y = z = will be called O. The geodesic distance 
between O and {x,y,z) is Rw. 

The part Pq of the pressure can be evaluated at O, i.e. the electric field in (3.2) is 
the one at O. Pq is that part of Pgeif which is created by the pseudocharges located at 
a geodesic distance from O smaller than RipQ. The regularization prescription is that 
there is no particle in a thin slab \x\ < e. The electric field E{w)t created at O by a 



15 



pseudocharge at {x, y, z) is given by (5.2) where = w and t = [—x/r, —yl'^^ —z/r). 
Thus, with (B.2) taken into account, the analog of (3.10) is 

:\x\<ro Jo cosw x"^ + p2 

where Tq = i^sin'^o- An expansion in powers oir/R gives 



Po^-—hm / —— — ^E^(w) (B.3) 

StT £-»0 y£<|x|<r-o ^0 COS W X^ + p^ 



^ = rT^[l + 0(rVi^^)] (B.4) 

COSW [x^ + p^)^ 

When the expansion (B.4) is used in (B.3), the leading term of (B.3) is (3.10) and 
the next term (which gives a convergent integral for which the e regularization is 
superfluous) is O{ro/B?). Using ro = i^sin'^o, one does find 

,2 



i.e.(5.8). 
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